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We consider a variational approach to the finite temperature Yang-Mills theory in the Coulomb 
gauge. The partition function is computed in the ensemble of glueballs and quasi-gluons which 
emerge as eigenstates of the Coulomb gauge Hamiltonian. We compute the energy density and 
pressure and compare with results of lattice computations for both SU (2) and SU(3). The emergence 
of a phase transition is discussed. 
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I. INTRODUCTION 

In recent years there has been an expansion in stud- 
ies aimed at the determination of the patterns of the 
QCD phase transitions [HQ. These arc crucial for an 
understanding of the mechanism of confinement and dy- 
namical chiral symmetry breaking. At high temperature 
and/or density, due to the asymptotic freedom, it is ex- 
pected that the weak interaction between quarks and glu- 
ons determine the properties of the quark-gluon plasma 
(QGP) 043- Lattice simulations at finite temperature 
are a good tool to investigate these phase transitions @- 
[l2j , while phenomenological models also enable studies of 
high density regime where the restoration of chiral sym- 
metry is expected [l3l - fl7j . 

The present paper investigates the thermal properties 
of a phenomenological model motivated by the canon- 
ical approach to QCD in the physical, Coulomb gauge 
quantization. We compute the partition function in the 
ensemble of glueballs and quasi-gluons. There are nu- 
merous studies of QCD thermal properties in covariant 
gauges that include for example Dyson-Schwinger based 
models, (l8l - [2l| or approaches based on the renormaliza- 
tion group flow 1221 . or direct models of the equation of 
state (EOS) |23h27| |. The few approaches that exist in 
physical gauges are rather loosely related to the underly- 
ing QCD interactions p^ - l32] . Recently, there has been 
also an attempt to explore the dynamical breaking of 
chiral symmetry in a self-consistent calculation at finite 
density [33|. 

The advantages of physical gauges for phenomenology 
and for developing physical intuition are clear, and we 
summarize them here. The degrees of freedom of the pure 
Yang-Mills (YM) theory are transverse gluons, and ther- 
mal excitations connect color-singlet states of arbitrary 
number of gluons. Transverse gluons are expected to be 
effective only at high temperatures, while at low temper- 
atures it is more effective to compute the partition func- 
tion in terms of the ground state glueballs [U, Hf| . The 
underlying interactions in Coulomb gauge are dominated 
by the instantaneous Coulomb potential acting between 



color charges. In the non-Abelian theory, the potential 
not only couples charges but it also depends on the gluon 
distribution of the state in which it is calculated. At zero 
temperature, in the vacuum state this distribution is such 
that the Coulomb potential becomes confining, i.e., pro- 
portional to the distance R between the external color 
charges, V(R) = acR HUHzl- Using various approxi- 
mate, variational models for the ground state YM wave 
functional, it has been possible to obtain a potential that 
is confining [38[ or almost confining, i.e., V(R) — > i? 1_e 
with e w 0(10%) [H-[4l|]. The Coulomb string tension a c 
is larger than the string tension computed from the tem- 
poral Wilson loop. This is because the Coulomb poten- 
tial represents the energy of a static quark-antiquark pair 
submersed in the QCD vacuum, while the Wilson loop 
measures the energy of the exact quark-antiquark state in 
which the gluon distribution is squeezed by closed vortex 
lines. Since the Coulomb potential is an instantaneous 
observable, one might expect that it remains confining 
even in the high-temperature limit [37|: At high tem- 
peratures the integration over transverse fields becomes 
even less restricted than in the vacuum, and, according 
to the Gribov-Zwanzigcr confinement scenario 

0,13, 

Coulomb confinement originates from large field configu- 
rations near the Gribov horizon. 

Recently the variational approaches to Yang-Mills the- 
ory in Coulomb gauge, [HI, Ho, 5l| have been extended to 
full QCD [44[ an d t° finite temperatures, 0, [|| assuming 
a quasi-particlc picture for the gluons. In the present pa- 
per we study the thermodynamic properties not only of a 
system of quasi-gluons but also include glueballs, which 
are the physical constituents of the Yang-Mills ensem- 
ble in the confining phase. In the following, we investi- 
gate the finite-temperature properties of Coulomb gauge 
Yang-Mills theory with focus on the aspects of thermo- 
dynamical properties and their behavior around the crit- 
ical temperature of the phase transition. In particular, 
we compute the energy density and pressure in the en- 
sembles of glueballs and quasi-gluons and compare with 
SU(2) and SU(3) lattice results. In Sec. II we present 
the general setting for the finite-temperature, canonical 
Coulomb gauge problem as well as general properties of 
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the thermal average when the glueballs basis is used. In 
Sec. Ill we present the numerical results. The summary 
and outlook are given in Sec. IV. Finally, details of the 
important expressions are presented in the Appendices. 



II. HAMILTONIAN APPROACH AT FINITE 
TEMPERATURES 

The Coulomb gauge Yang-Mills Hamiltonian is ob- 
tained after gauge fixing and elimination of the Gauss's 
law constraint on the longitudinal component of the elec- 
tric field, 



H YM = i J [3}x (J- 1 [A] UJ[A] n + B 2 ) + H c 



(1) 



H 



K 



He 



H c 



[3M% J- 1 [A] P a (x) J[A] Ff(x, y) P \y). 



(2) 

Here II a (a:) = — \8/8A a (x) is the canonical momentum 
(electric field) operator, and 

J [A] = Det(-DV) (3) 

is the Faddeev-Popov determinant with 

D = V+gA, A ab = f c A c , (f c ) ab = f acb (4) 

being the covariant derivative in the adjoint representa- 
tion. Furthermore, 

p a (x) = -f abc A b ■ IF (5) 

is the color charge density of the gluons and 

Ff(x,y) = (X.OK-2JV)- 1 (-V^i-DVy^b) (6) 

is the gluon field dependent Coulomb kernel. The vac- 
uum expectation value of this kernel plays the role of an 
instantaneous potential between color charges. At zero- 
temperature and at large distances (Fjj"(x, y)), is well ap- 
proximated by a linear rising potential. The gauge fixed 
Hamiltonian in Eq. ([TJ is highly non-local due to Coulomb 
kernel Fa{x, y), Eq.©, and the Faddeev-Popov determi- 
nant J[A], Eq.Q. In addition, the latter also occurs in 
the functional integration measure of the scalar product 
of the Coulomb gauge wave functionals 



(Vi|C#2) = / DAJ[A]il>l[A]Oih[A]. 



In Ref. [40( the Yang-Mills Schrodingcr equation was 
solved by the variational principle using the following 
ansatz for the vacuum wave functional 



(A\0) = 



1 



(A\0), 



(A\Q)=Afexp[-- t 



A(-k)uj(k)A(k) 



(8) 



The pre-exponential factor removes the Faddev-Popov 
determinant from the scalar product Eq.©. The kernel 
ui(k) was determined by minimizing the vacuum energy 
(H Y m), which yields an w(fc) that can be well fitted by 
Gribov's formula 



M 4 
V 2 



(9) 



and which is in satisfactory agreement with the lattice 
data (HI, for M « 880 MeV. 

The present paper is devoted to study Yang-Mills the- 
ory at finite temperatures, which is defined by the density 
operator 

V = Z- 1 exp(-pH YM ) 1 (10) 

where /? = 1/T is the inverse temperature and 

Z = Tre- /3ff ™ (11) 

is the partition function. Exact calculation of the trace 
in the thermal averages 



(O) = Tr(OD) 



(12) 



is not possible for the YM theory and the way we pro- 
ceed is, following ref. [l| to replace V by a variational 
ansatz. This is achieved by first defining a suitable basis 
in the gluonic Fock space. It is chosen as follows. In the 
standard fashion we Fourier decompose the gauge field 
in terms of creation and annihilation operators 



4 Q (x) 



n«(x) = -i 



a?( q ) + a?(-q)}e-^ 



(2tt)3 y ^J(qj 
(2tt) 3 



K(q)-af(-q)]e-^ 

(13) 



where a b (q) = J2x e »(9> ^) a ('Z: ^; ^ are tr± e helic- 

ity and color indices, respectively). Choosing here u>(q) 
to be the kernel of the vacuum wave functional © the 
operators af(q) annihilate this state, i.e., 



<(fc)|0)=0. 



(14) 



Then a complete basis in the gluonic Fock space is given 
by 



{|n}} = {|6>, af(k)\Q), af(k)ay(q)\b), ...}. (15) 

As shown in Re.[l|, however this quasi-gluon basis, even 
when restricted to color singlet states is not ideal for stud- 
ies of thermal properties of the YM plasma. Because of 
confinement, energies of isolated gluons are large (infi- 
nite) and subtle cancelation of infrared divergencies have 
to occur in color singlet states containing multiple gluons 
with low relative momenta. For the same reason gluons 
in these low momentum states are expected to strongly 
bind into color singlet, aka glueball states. Thus instead 
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of working directly with the basis of Eq. (fT5"j) we choose 
to work with a basis of glucballs, which are constructed 
using creation operators defined by, 

Gt («) = yE r ( 1 ' 2 ) flt ( 1 ) at ( 2 )' ( 16 ) 

(1,2) 

where V = (2tt) 3 S(0) is the volume factor. Here the 
summation extends over single gluon helicitics (Ai^) and 
color (ci^) and it also includes integration over individ- 
ual gluon momenta, (pi 2). When acting on a vacuum 
state creates a glueball state with quantum numbers 
a = (P, J PC ) where P and J PC are the total momentum, 
and the total spin, parity and charge conjugation of the 
glueball, respectively. The cut-off on relative momentum 
is implicit in the glueball wave function, \& a . The latter 
can be written explicitly as 



^(1,2) = ^ 



A2;ClC2 



(Pl,P2) 



(2nf8 i (P-p 1 - P2 ) 



and it is normalized by 



V2 



A1A2 



[dpidp 2 ]p4 , A 1 A 2 (Pl^2)| 2 



(17) 



(18) 



with the measure defined as, 



[dp 1 dp 2 ]p = 



d 3 pi d 3 p 2 
(2tt) 3 (2tt) 3 



(2Tr) 3 5(P- Pl -p 2 ). (19) 



Bose symmetry implies that the wave function is sym- 
metric under exchange of the quantum numbers of the 

two gluons, *£ lAa (Pi,p 2 ) = ^x 2 \ l (P2^Pi)- In tcrms 
of this single glueball operator, multiple glueball states 



are given by 



|71 Ql Tl a . 2 



n 



-|0>. 



(20) 



In the following we will ignore the the Faddeev-Popov 
(FP) determinant. As will be shown below, the expres- 
sion we obtain from variational principle is closely related 
to that for the spectrum of Hym and at the end given by 
the eigenvalues of the Hamiltonian. Thus the FP contri- 
butions to the formulas for the free energy can in princi- 
ple be restored by comparing with those for Hym ■ 

We introduce a variational ansatz for the thermal den- 
sity operator by replacing Hym in Eq. (fTU)) by the follow- 
ing single particle operator 



(27T) 



(21) 



The optimal value of free energy is obtained by taking 
a variation with respect to fi, which obviously has the 
meaning of the gluon energy. 



A. The partition function 

Computation of the partition function (jllj) with Hym 
replaced by h (|2~Tj) in the glueball basis ((20)) 



Z= (n ai n a2 ..\e-^^ n ^\n ai n a2 ..), (22) 



is straightforward and yields 
d 3 P 



exp 



V 



n a (P)) 



where 



n a (P) 



-0E a 



1 



(23) 



(24) 



is the glueball thermal occupation number. Here the ef- 
fective glueball energy, E a given by 



-PE a (P) 



= J2 fidpid P 2}p\n 1 xM>P2)\ 2 e- 0nir - 



where 



fii, 2 = n(pi) + n(p2). 



(25) 



(26) 



Since the density operator is defined in the gluon basis 
but the thermal average is evaluated in the basis of gluc- 
balls the effective Boltzman factor, exp(— j3E a ) is deter- 
mined by averaging gluon thermal distribution over the 
glueball wave function. 



B. The internal energy 

Using the Fourier decomposition (fTB"]) we can calculate 
thermal expectation values of the Yang-Mills Hamilto- 
nian Eq.(TT]). After normal ordering, the Hamiltonian 
contains the vacuum contribution and one- and two- 
body gluon operators. Those are given explicitly in Ap- 
pendix [X] Since vacuum contribution is temperature in- 
dependent it can be removed by defining the free energy 
with respect to that of the vacuum. The final expression 
for the thermal average of the Hamiltonian is then given 
by 



(H 



YM 



V 



E 

JPC 



d 3 P 

—lE a (P)+B a (P)+C a (P)][l+n a (P)} 



(27) 

where the three terms represent contributions from: the 
the one-body operators describing single gluon ener- 
gies averaged over the glueball state, (£), the two-body 
magnetic contribution from the four-gluon vertex (B), 
and the two-body Coulomb interaction, (C), respectively. 
The explicit formulas for the three terms are given in the 
Appendix [5] 
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C. The free energy 

The glueball wave function is the solution of the Hamil- 
tonian bound state problem and as discussed previously 
defines the basis over which thermal averages are com- 
puted. The variational estimate for the free energy, F 



T = (H- 



YM) 



TS 



(H 



YAI 



InZ dlnZ 



(28) 



is in turn obtained by minimization with respect to the 
single gluon energy tt(k). Before proceeding, however, 
we note that boost invariance requires the density ma- 
trix to depend only on the total momentum of the two- 
gluon state. Thus the factor exp(— /3(Q(pi) + fi(p 2 ))) 
which appears in matrix elements should be replaced by 
cxp(— /3a/ + M^) where M a is a Lorentz scalar. From 
Eg. (1231) it then follows immediately that, 



E a (P) = VPj+M* (29) 
and the internal energy given by Eq. ([2"T]) becomes 
(H YM ) V f d 3 P 

JPC ■ 



(2tt)' 



■n a (P)[C(P)+B° a (P)+C° a (P)} 



(30) 

where the subscript indicates that the correspond- 
ing quantities are to be evaluated at j3 = (of. 
Eqs. (|B1 ^ .([B"T 1) .([B"8 |) ). Finally the entropy reduces to Q 



S 
V 



E 

JPC 



d 3 p 

(2tt) 3 



[ln(l + n a (P))+pE a (P)n a (P)}. (31) 



Minimization of F is now performed with respect to M a 
and yields the following relation, 



£ a (P) + B a (P)+C a (P) = E a (P) 



(32) 



which is immediately recognized as the zero-tempcraturc 
eigenvalue equation for the Hamiltonian Hym projected 
onto the glueball sector. It thus follows that M a is the 
glueball masse obtained from diagonalizing the Hamilto- 
nian. With eq. ([32]) we find from eqs. ([30]) and §31$ for 
the minimum of the free energy 



F 

V 



-y 

JPC 



d 3 p 

(2tt) 3 



In 1 



(33) 



D. Scattering states 



becomes relevant. At this point we qualitatively describe 
the spectrum in terms of single quasi-gluons, which are 
no longer confined but screened. Calculation of the ex- 
pectation value of Hym in the basis f| 15[) of quasi-free 
gluons was done in and gives, 



(H 



YM) 



V 



2{N 2 C - 1) 



d 3 q 
(2tt) 3 



n(q)e(q) 



d 3 q d 3 p 
(2tt) 3 (2tt) 3 



[Hp, q) + c (p, i)] n (p) I 1 + n (i)} ■ 

(34) 

Here the three terms, given explicitly in the Appendix [B] 
represent contributions from the single gluon energy, the 
four-gluon magnetic term and the Coulomb interaction, 
respectively. The free energyis obtained from Ea. ([28]) 
with gluon entropy given by [l| 

| = 2(A£ - 1) / (0 Ml + »(«)] + m<i)<q)] ■ 

(35) 

Minimizing the free energy with respect to the density 
matrix SF/5fl(k) = results in the following expression 
for the effective gluon energy il(k) 



n(k) = e(q) 



d 3 p 

(X) 3 



[b(q,p) + c(q,p)}[l + 2n(q)} (36) 



which, when substituted into the expression for free en- 
ergy, yields, 



F 
V 



2(JV£-1) f d 3 q 







2{Nl 



(2^) 3 
d 3 q d 3 p 
(2tt) 3 (2tt) 3 



ln[l + n(q)} 



4n(p)n(q)[b(p, q) 



c(p, q)]. 
(37) 



Here the first term represents contribution from the free 
gas of quasi gluons and the second one is the one-loop 
correction due to residual interactions. We emphasize 
that in order for the two-component glueball plus gluon 
model to be valid, the interaction in the Coulomb term in 
Eq. Q37[) (implicit in the term proportional to c) has to be 
screened. That is, c(p,p) is assumed to be free from the 
infrared singularity at p = q normally associated with 
confinement. Scattering corrections, i.e. the one- loop 
term are then expected to be week and we ignore them 
in the numerical studies. 



If the Coulomb interaction, F^ b (x,y) (c.f. Eq.([T])) is 
replaced by its vacuum expectation value, which is ex- 
pected to be strictly confining, [36j], quasi-gluon bound 
states saturate the spectrum of Hym- However, since the 
Coulomb kernel couples gluon Fock sectors with arbitrary 
large number of gluons, reduction of the full Hamilto- 
nian to the two-gluon subspace must break down above 
energies where excitations of multiple quasi-gluon states 



III. NUMERICAL RESULTS 

We present a thermodynamical study for the energy 
density and pressure of the SU(2) and SU(3) gauge the- 
ories separately for the quasi-gluon and glueball ensem- 
bles. In the glueball ensemble the degeneracy factor de- 
termines the high temperature limiting value of the en- 
ergy density but the behavior of the transition is not 
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well know as well as the location of the critical temper- 
ature. In the thermodynamical quasi-gluon study, the 
phase transition, the critical temperature Tc, the behav- 
ior below and above the transition and the high temper- 
ature limit are compared with SU(2) and SU(3)-lattice 
results (23T - [26l | . when two dispertion relations f2(fc) are 
used. 



In the case of glueballs we need to know the expected 
degeneracy. Since explicit digitalization of the Coulomb 
gauge YM Hamiltonian [34|, |35| reproduces lattice glue- 
ball spectrum [49| we use the latter to determine the 
number of states. In particular we will consider glueballs 
up to 2.5 GeV, H^fl], i.e. with J PC = 0+-,(T+,2++. 
The numerical results for energy density and pressure, are 
shown in Figs. [TJ[5J At high-temperatures T >>,M a 



A. Glueball energy density and pressure 

Energy density and pressure are computed from the 
free energy differentiating In Z with respect to T and V, 



T 2 d\nZ 



V dT 



p = T 



d\nZ 
dV 



(38) 



(39) 



with Z = exp(— /3J-). The energy density in the pure 
gauge theory has been found to rise rapidly at Tc 
and approach the high temperature ideal gas (Stcfan- 
Boltzmann) limit from below (24) . In the high tempera- 
ture limit one finds in leading order perturbation theory 
for SU(N) H3 
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FIG. 2. Same as in Fig. fj] for the combination of energy 
density and pressure (e — 3p)/T 4 . 
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(42) 



where M = is the degeneracy factor. Inspecting the 
numerical results it is evident that because of the high 
degeneracy, glueballs contribute to much to thermody- 
namical quantities as compared with QCD expectations 
from lattice simulations. It implies that glueballs must 
evaporate below the critical temperature, and a Hamil- 
tonian model based on confined potential without mix- 
ing with open-channels becomes inadequate at fairly low 
temperatures. 



FIG. 1. Energy density versus temperature in the case of 
glueballs for J PC = 0+-,(T + ,2 ++ . 

In the following we denote energy density and pressure 
for the ensemble of glueballs by e\ , p\ and for the en- 
semble of gluons by 62 , P2, respectively. In particular, 
glueball energy density and pressure are given by 




B. Gluon energy density and pressure 

The gluon energy density and pressure are given by 

„ = 2[ivS-i]r/|L ln ( 1 + _J r _) 

which, assuming Q,(q) — > q for large moments, in the high 
temperature limit reduce to 



6 



2(N C - 1) T 4 
2(7V2 - 1) T 4 1 



£2 | _ 7T 

T->oo — 



2 

30' 



7T 

90' 



(44) 



which, of course, yields the correct Stefan-Boltzmann 
limit given by eq. (|40[) with a s = 0. 
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FIG. 3. Energy density versus temperature for Nc = 2 gluon 
ensemble with Gribov dispersion relation. 
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FIG. 4. Same as in Fig. [3] for the combination of energy 
density and pressure, (e — 3p)/T 4 . 

For the gluon dispersion relation we use the Gribov 
formula 



m„ 



(45) 



and choose the Gribov mass m g in the range from zero 
(perturbative gluons) to 880AfeV '. The latter value is 
found on the lattice [H[. In Fig s - LH EH an d Figs. [5] we 



summarize the results for Nc = 2 and Nc = 3, re- 
spectively. Of course, in this calculation the transition 
temperature is set by the Gribov mass m„. For the lat- 



tice value 



880Mey the quasi-gluon ensemble re- 



produces the lattice energy density reasonably well up to 
the transition temperature but substantial deviations oc- 
cure above the phase transition. This should come with 
no surprise. In a self-consistent treatment of the finite 
temperature quasi-gluon ensemble in the variational ap- 
proach in Coulomb gauge [|J one finds that at the decon- 
finement phase transition the gluon dispersion relation 
switches from the Gribov formula (|45[) in the confining 
phase to the massive dispersion relation. 



Q(k) = \/k 2 + to 2 



(46) 
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where the gluon mass m is temperature dependent and 
growths linearly in T for large T. Its minimal value 
can be as low as 200MeV. We have also computed en- 
ergy and pressure for the massive dispersion relation (|46[) , 
with m in a range from zero (perturbative) to 880 MeV. 
As expected, the Gribov formula reproduces the critical 
temperature and the overall shape more accurately than 
the massive dispersion relation. 



IV. SUMMARY AND OUTLOOK 

We studied the Coulomb gauge Yang-Mills theory at 
finite temperatures using a variational approach. The 
partition function has been computed in the ensemble of 
glueballs and quasi-gluons. Working with both ensembles 
we present the possibility of a phase transition since glu- 
ons with low relative momenta are expected to strongly 
bind into color singlet, glueball states. The thermody- 
namical limits for the energy density and pressure are 
different in each ensemble. This is expected since the par- 
tition function of the glueballs depends on the states J PC 
included and their degeneracy, while the partition func- 
tion of the (quasi-)gluons depends on the number of col- 
ors Nc ■ In the present work we have considered the glue- 
ball states up to M G = 2.5 GcV, i.e. J PC = 0~, 0",2 
and showed that the thermodynamical limit is rapidly 
overshoot, indicating the possibility that the glueballs 
may evaporated at some finite temperature. Further- 



more, the more glueballs one adds the higher the degen- 
eracy factor and the thermodynamical limit will be over- 
shoot faster, indicating that the glueballs evaporate at 
quiet low (below Tq) temperatures. A more realistic de- 
scription of the deconfinement phase transition would as- 
sume a two-component picture, in which glueballs coexist 
with gluons. Well below the deconfinement phase tran- 
sition the finite temperature Yang-Mills ensemble would 
dominantly consist of glueballs, which dissociate at the 
deconfinement phase transition into pairs of gluons. 
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Appendix A: YM Hamiltonian Contributions 

The Hamiltonian thermal average is computed in the 
ensemble of glueballs and quasi-gluons which immedi- 
ately rules out contributions from an odd number of par- 
ticle creation/annihilation operators. The relevant terms 
of the YM Hamiltonian are 



Hk \ I (0 Mp ? Vu(qM q i)(2nf6(q + q')[a?(q)a? (-q') + a? (-q)a?(q')} , (Al) 

a.i 

A J (2nr (2wr ^ vM<zM<7') 

g 2 f d 3 q 1 d 3 q 2 d 3 q 3 d 3 q 4 ^ ^ (2 7 r) 3 (5(q 1 + q 2 + q 3 + q 4 ) 

+ 32 J (2n) 3 (2n) 3 (2,) 3 (2*)S 2^2^ »* ilm ^iM^-teM^) 

'a% qi )al(q 2 )af(- q3 )a^(- qi ) + afajctf \-q 2 )af(q 3 )a%(-q 4 ) + afajag (-q 2 )af (-q 3 )a e m (q 4 ) 

+ a b ^- qi )al{ q2 )ai{q 3 )a^{- qi ) + af (- qi )a c k (q 2 )af (-q 3 )a e m (q 4 ) + of (-qja* (-q a )af (q 3 )a e m (q 4 ) 

(A2) 

„ ,9 2 f d 3 qi d 3 q 2 d 3 q 3 d 3 q 4 ^ ^ bc ade (0 s3 r/„ +„\J?f „ „\ MggM^j 

a»(g^(g 2 )af(-g 3 )af(-^^ 
+^(-q 1 )a c M^( q3 )af(- q4 ) af(- Ql )^(g 2 )af (-qjafa) + a* (- Ql )af (-q 2 )<(g 3 )a J e (q 4 )" 

(A3) 



To compute the thermal averages we need to write Computation of thermal averages in the gluon basis (|15[) 
products of particle operators in normal ordered form. 
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was given in [l| while in the case of the glueball ensemble the relevant matrix elements to compute are 
I 



<^>^</^(V[^-l]^+E^«r(«K(ff))> (A4) 



(Hb) -+ {\ I ^ J2 ^ilm^(Vm - l]t km {q) + (<l)<&(<l) + (k O m)]) 

3 2 /" d 3 g d 3 q' ^ f^,, AT 2 ,,,, tji{q')t km {q) - t k i(q')t jm (q) 



kI<^U*^ eijkeiim { v[N2c - i]Nc 

ijklm 



u(q)w{q') 



J2 e abc e ade ^ *M (g ) ^c^ef ( g ) Q b ( g ) + (fc ^ £ 6 ^ c ) + ( m o e) + (j fl M« « e)] 

abode 

g 2 f d 3 qi d 3 q 2 d 3 q 3 d 3 q 4 ^ ^ , a 6c «d e 



' 7 (2n) 3 (2n) 3 (2^ (2,) 3 2^ h_ ' 



32 J (2n) 3 (2n) 3 (2n) 3 (2n) 3 ^ «* y/ufaMtoMtoMto) 

[(2n) 3 S 3 ( qi + q 3 q 2 q 4 )a$ (q 2 )a%(q 4 )a b A qi )at (q 3 ) + (3 O 4) + (1 O 2) + (1 o 2; 3 O 4)] ) (A5) 



■1 + ^ 
toq 



tij(q)tij(q') 



£ e abc e ^2^(g')F( g - g')^[^ d af (g)a|(g) + (bo c,d o e)] j 

abede / 

f d3qi d3q2 d3q3 d3qA v" v abc ade / ^ggMgjj 

8 7 (2tt)3 (2tt)3 (2tt)3 (2tt)3 ^j^/ ' V w(gi)w(g 3 ) 

[-(27r) 3 5 3 ( 9l + q 3 - q 2 - 9 4 )i ? (-93 + ^K^^teJaJtei W (<7 3 ) + (3h4)|(1h2)-(1h2;3o4)]) 

(A6) 



with further details given below. Appendix B: Hamiltonian Thermal Average 
I 



1. (Hym) in a basis of glueballs 

Computation of thermal averages of the Hamiltonian in the glueball ensemble involves expectation values of one- 
body and two-body operators. These are given below 



d 3 q ,, s h t 



f( q )a?(q)a%q)) 

d 3 P a d 3 q r , J . rj/j/1 <5 &c e -« fi «}+ fi «» 



Nq — 1 1 - e-^ B ° 



Jf c Ai ,A 2 ,Ai ' 

x [(2^) 3 «5(p 2 - p 2 )(27r) 3 % - p[)(27:) 3 S(q - pi)/(g)*£ liA| (pi^^^^^.a.Cp'i.^) + (pWp' 2 )] 
d 3p jhc e -/9(n(pi)+n(pi)) 



V E E I^wm p 



iV£ - 1 1 - e -^o 



J pc Ai,A 2 ,Ai ' 

xIMW,,^ A^.Pl)*U{Pi.Pi) + W ( B1 ) 

I 

where (q, q) = H 1 ^ ^ (q) — e*(q, Xi)ej(q, X 2 ) is the In a similar way the two-body operator thermal average 

one-body vertex factor and 1/1 — e^ l3Ea = [1 + n a (P)]. 
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is given by 



V 



d 3 qi d 3 q 2 d 3 q 3 rf 3 q 4 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 

d 3 P Q 



/(<7i,'72,'73,'?4)(27r) 3 (5(g 3 + g 4 - gi - q^a* (gi>f ^^o^^K 4 ^)) 



(27T) J 



[dpi dp 2 ] p [dpi dp 2 ] p X X 



-0(O(Pl)+ n &>2)) JC3C4JC1C2 

1 - e-/ 3 ^ iV^ - 1 

7 p c A1A2A3A4 G 



X {/bl,P2,pi,^)*? lA2 (pi,P2)<T 2 A 3 A 4 (Pl.P2,pi,^)^ 3A4 (pi,^) + (Pi ^P' 2 )} 

where the vertex factor for the two-body operator is given by 

HZ^XsxSP^'P^P'nPz) = e **(Pi> A i) e j(P 2 ' A 2) £ r(pi, A 3 )e s (p 2 , A 4 ). 

I 



(B2) 
(B3) 



Calculation of the thermal average of the Hamiltonian This finally leads to the expression in Eq. (|27p with 
is further simplified by the following relations involving 
vertex factors 



E^(9 , )^fA 2 (9)=2 ( 5A 1 A 2 

jk 



Ai .A2 



£ »(P) = E / [dpid P 2]p\n u xM^ 



E^(9')<a 2 (9) = E 

jk j 



2 H{ J iX2 {q,q) 



OAiA, = " 



Ai A 2 



2 ^ 2 2 

E^1a 3 A.(^-9'>9>^-9.9') 

ij 

ij 

^2 H \iiU 3 xA Pa ~ q',q',q,P a - q) « <5 AiA2 6a 3 a4 



x [e(p 1 )+e(p 2 )]e- /3 ( n( f 1 ) +n (f 2 » (B5) 

with the single gluon energies, e(p) given by a sum of 
kinetic and self energy terms, e(p) = oj(p)/2+p 2 /2w(p) + 
£ B (p) + £ c (p) 



5 2 iV c f d 3 q 3 



8 7 (2tt) 3 w(p)w(g) 



A2 A4 



The explicit form of the two-body magnetic (B) con- 
tribution from the four-gluon vertex and the two-body 
(B4) Coulomb interaction, (C), respectively are 



J 



B a (P) 



9 2 N C 
8 



E 



[dpidp 2 ]p[dpidp 2 ]p 
\/cj(PiMP2MPiMp 2 ) 



-/3(f2(pi)+f2(P2)) 



*A 1 A>l.P2)*V 1 A 2 (Pl.P2)(^ 1 A 2 5AiA 2 - ^V^AaAj) + fri Pa) 

I 

I 



(B7) 



e a (P) 



g 2 N c 



A, 



E / [dpidp 2 } P [dp[dp' 2 } P e-^ n ^ +n ^ [F( Pl -p' 1 )n i x 2 iPuP^n i xM^2) 



/ /cj(pi)o;(p 2 ) uj{p 1 )uj(p2) \ , , , , 



(B8) 



2. (Hym) in a basis of quasi-gluons result 



The calculation of the thermal average (Hym) in the V J (27r) 

basis (|T5j) of quasi-gluons was carried out in [l|with the /■ 



2(^0 - 1) / 7^kn{q)e{q) 

b(p, q) + c(p, q)]n(p) [1 + n(q) 



(2tt) 3 (2tt) 3 1 

(B9) 
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with 



b(p, q) 



u{q) q 2 



Mi) 



g 2 N c 3-p-q 
8 oj(p)uj(q) 



(BIO) 



(Bll) 



c(p, q) 



— [l + (p. 9 ) 



(B12) 
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